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Abstract 

We show that the third-order difference equations proposed by Hirota, Kimura and Yahagi are 
generated by a pair of second-order difference equations. In some cases, the pair of the second-order 
equations are equivalent to the Quispel-Robert-Thomson(QRT) system, but in the other cases, they 
are irrelevant to the QRT system. We also discuss an ultradiscretization of the equations. 



PACS numbers: 02.301k, 05.45.-a 



I. INTRODUCTION 



Discrete integrable systems have attracted much attention and a lot of studies have 
been done by various points of view, such as integrability criteria (singularity confmemen ; 
property l|, algebraic entropyj^J), geometric or algebraic description of the equations^, I4, 
110,0, E| and so on. 

Especially, second-order integrable difference equations including QRT system P. Iic| and 
discrete Painleve equations which are regarded as non-autonomous variations of QRT 
system, have been extensively studied, and a number of significant properties have been 
obtained. 

For example, a symmetric version of QRT system is defined by the following form 
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with arbitrary symmetric 3x3 matrices A and B. This equation has a conserved quantity 
h(x n -i,x n ) defined by A and B. Moreover, a general solution is described by an elliptic 
function. 

However very few results have been obtained for third-order integrable difference equa- 
tions. Research on such equations are important in order to reveal integrable structures of 
general discrete integrable systems. 

In this paper, we investigate third-order integrable difference equations proposed by Hi- 
rota, Kimura and Yahagi 12( and show that they are generated by a pair of second-order 
integrable difference equations. Moreover we also discuss their ultradiscretization. 

Hirota, Kimura and Yahagi have investigated third order difference equations of the form 
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and have found that nine equations, 
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are integrable in the sense that they have two independent conserved quantities. A re- 
markable property of these equations is that their trajectory of a solution in 3- dimensional 

phase space looks like a composition of two separate curves. Fig. ^ is an example of such 
trajectories in 3D phase space which is generated by 

y n+ 2y n -i = a + y n + y n+1 , (4) 

where the equation is obtained through a variable transformation y n = — — , a = 

b\x n &i 

from Eq. (1Y6I) . Moreover, it is an important fact that odd step points belong to the one 



curve and even step points belong to the other. 

This fact strongly suggests that a combination of lower dimensional integrable equations 
determines the integrability of the third-order difference equation. We show this is true for 
all nine equations in the following section. 
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FIG. 1: A trajectory of a solution to Eq. Q for a = 2.0, yo = 1.0, yi = 2.0, 1/2 = 1.5 

II. PAIR OF SECOND-ORDER INTEGRABLE EQUATIONS GENERATING A 
THIRD-ORDER EQUATION 

A. Y6 



If we take a backward difference of Eq. P|) 

y n +2Vn-i = a + y n + y n+ i, 

we obtain 

A„(y n+2 Z/ n .-i -a-y n - y n+1 ) = 0, (5) 
where A„ is a difference operator defined by A n /„ = f n — f n -i- Eq. (0) can be written as 

(1 +y n +2)(l +l/n) = (1 +l/n)(l + l/n-2) 
Vn+l Vn-l 

This formula means that there are constants which depend on the initial values and on a 
parity of n. Hence, we obtain 
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where # n = y 2n , K = y 2n +i, c = , c x = . From Eqs. (J4J) 
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and (J7J), we obtain a pair of second-order difference equations 
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where Eq. (JHJ) is a equation for even steps and Eq. (JHJ) is a equation for odd steps respectively. 
Eq. (jHJ) and Eq. Q can be written in the QRT form 
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Consequently, conserved quantities of Eqs. (jHJ) and © are given as 

h =(glgl+i + (2 + c )(g n + g n+1 )g n g n+1 + (1 + c )(g 2 n + g 2 n+l ) 
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h =[Kh n +i + (2 + ci)(/i n + h n+l )h n h n+1 + (1 + Cl )(/i 2 + u n+1 ) 

+ (2 + 2ci + a Cl + c 2 )(/i n + h n+1 ) + (1 + d)(l + a Cl ))/(/i„/i n+1 ). (14) 

Hence invariant curves of Eqs. (JHJ) and (JHJ) are given by the above equations. These curves 
determine the structure of the trajectory of a solution to Eq. (|Y6J) in 3D phase space shown 
in Fig. ^ This is the simplified integrability structure of Eq. (JY6|) and we show below a 
similar structure exists in the other eight equations. 
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From Eq. (jYT 
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Xi 
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^2 

From Eqs. (jYl|) and (JT7J), we obtain a pair of QRT systems 

Gi{g n ) - 9n-iG 2 {g n ) 



9n+l 



) _ 9n-iG 3 (g n ) 

H\{h n ) — fe n _iiJ 2 (fe n ) 
H^ihn) — h n _ 1 H 3 (h n ) ' 



where 



A(c) 



b{ - 2a 3 b 3 + 







a\ — a\bi 



\ 



2a\ — a b 3 + 2a 3 c — 



dib 3 C 

~bV 



B{c) 



h h 

61 -c 




3 - a l b l 

\ 



2a\a 3 — CL^bi + a±c 



2a\a 3 — a^bi + a,\C 



a\ + a c J 




7 



where 



A(c) 



«22 
023 



/ 



(a — ai — c)(ai — 1) — (a — a x — c)(ai — l)a 2 
a 22 a 23 

,2 



y-(a - ai - c)(ai - l)a 2 a 32 a 33 y 

2(a — — l)a 2 + c(a — 2ai + — c), 

2(a — a-i)(ai — l)a^ + cai(— ao + 2a oi — 2a 2 + — aic) 



a 3 2 — «2 3 , 



a 33 = — (a — a\)(ai — + cai(—a ai + 2a a 2 — a\ — Oqc) 



(29) 



' 1 ai 0^ 



5(c) 



a! a 2 — c 




(30) 















( H 1 {x) y 




f x 2 ^ 




^x 2 ^ 


G 2 (x) 


= A(cb) 




xB(c ) 






H 2 (x) 


= A( Cl ) 


X 


xB(ci) 


X 


G 3 (a:) y 




w 




lv 








v) 




Iv 



(31) 



D. Y5 



From Eq. (jY5 
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Note that we don't take a backward difference but a sum of Eq. (jY5|) here. Eq. 
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From this equation, we obtain 
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E. Y7 



Introducing a variable transformation x n = to Eq. ()Y7|) 
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F. Y8 



Introducing a dependent variable transformation x n = iy n , a[ = ia\ to Eq. (|Y8|) 

rp rp ■ ^ t ( Y R\ 
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Introducing a variable transformation x n = ^ n ~ l to Eq. (JY9 
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H. Y2 
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where g n = x 2n , K = x 2n +i and 

C = (1 + X )(l + Xi)(l + X 2 )(a(l + X + Xi+ X 2 ) + XqXi + X\X 2 + X 2 X + X XiX 2 ), 

X XiX 2 

(52) 

b = {(a - l)x x 2 (l + xif 

XqX~lX2 

+ (1 + xi + x xi + x 1 x 2 )(a(l + x + x\ + x 2 ) + xoXi + x\x 2 + x 2 x + x xix 2 ), 

(53) 

bi = — - — {(a - l)xix 3 (l + x 2 ) 2 
x x x 2 x 3 

+ (1 + x 2 + x\x 2 + x 2 x 3 )(a(l + xi + x 2 + x 3 ) + x\x 2 + x 2 x 3 + x 3 x x + xix 2 x 3 ). 

(54) 

Eq. (J5"T| does not belong to the QRT system as it stands, though there is still a possibility 
that it could be reduced to the system through a change of variable. 



I. Y3 



About Eq. (|Y3|) . we can only show numerical results. For various parameter do, a 3 and 
initial values xo ~ x 2 , it is generated by a pair of second-order equations on even steps (g n ) 
and odd steps (h n ). In any case, both equations follows 



_ 0<ij<3 , _ 0<i,j<3 ( v 

9n+i — x : r, n>n+i — x : r> l^J 

0<i,j<3 0<ij'<3 

where ~ cfjj are constant obtained from parameters and initial values numerically . This 
fact strongly suggests that Eq. (jY3|) is also derived from a pair of second-order equations. 



III. ULTRADISCRETIZATION 



In this section, we consider an ultradiscrete version of the third-order integrable 



equations 



13 



14j. Since ultradiscretization requires a positivity of parameters and depen- 



dent variables of the equations, (Y1),(Y2),(Y4),(Y6),(Y7) and (Y9) are ultradiscretizable. 
We show the procedure of ultradiscretization works well by Eq. (j3J) as an example. 
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If we use transformations y n = exp 



Y„ 



a = exp ( — ) for Eq. (|3J) 



y n+2 y n -i = a + y n + y n+u 



and take a limit e — > +0, then we have 



Y~ n+2 = max(A, Y~ n , Y~ n+1 ) - Y~ n _ 



(56) 



This is an ultradiscrete version of Eq. (Jl|. FigJ2I shows a trajectory of a solution to Eq. (|56j) 
in 3D phase space. 




FIG. 2: A trajectory of a solution to Eq. (fSfij) for A = 2.0, Y = 2.0, Yi = 2.01, F 2 = -5.1 



It follows from the result for Eq. (J3J) in the previous section that Eq. (J5T)|) is generated by 
a pair of ultradiscrete QRT system 



U n+ i = max(0, A + C , U n + max(0, C )) — U n -\ — max(0, U n ), 
V n+ i = max(0, A + Ci, V n + max(0, Ci)) - K-i - max(0, V n ), 

C = max(0, Y ) + max(0, Y 2 ) - Y u 

d = max(0, Yi) + max(0, F 3 ) - Y~ 2 , 



(57) 
(58) 
(59) 
(60) 
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and invariant curves become 

max(u n + £4 + i, max(U n , U n+X ) + max(0, C ), max(U n - U n+1 , U n+X - U n ) + max(0, C ), 
max(-U n , -U n+ i) + max(0, C , A + C , 2C ), 

-U n - U n+X + max(0, C ) + max(0, A + C )) = K , (61) 

and 

max(V n + K+i, max(y n , K+i) + max(0, C x ), max(K - Ki+i, V n+1 - K) + max(0, C x ), 
max(-14, -K+i) + max(0, C x , A + Ci, 2Ci), 

-V n - K+i + max(0, Ci) + max(0, A + C x fj = K x (62) 
FigEl shows invariant curves for Eq. (JoTj) determined by Eq. ()6ip. 




FIG. 3: An invariant curve of Eq. (|H7|) for A = 2.0, Yb = 2.0, Yi = 2.01, Y~ 2 = -5.1 

IV. CONCLUSION 

In this paper, we have shown the third-order integrable difference equations proposed 
by Hirota, Kimura and Yahagi are generated by a pair of second-order integrable differ- 
ence equations. In the case of Eqs. (|Y1|) and (|Y4JI ~ (IY9|) . second-order difference equations 



14 



are a special case of the QRT system. In the case of Eqs. (|Y2|) and (jY3|) . second-order 
equations may not be the QRT system. Furthermore we have shown the procedure of ultra- 
discretization works well for third-order equation, and derived second-order equations and 
invariants curves are also ultradiscretizable . 

Although the whole integrability structure of the general third-order equations is still 
unknown, our work would be one of keys to understand the structure. Generating our 
results, that is, investigating a connection between the general QRT system and the third- 
order equations is an important future problem. 

Acknowledgments 

The authors express our sincere thanks to Professor Ryogo Hirota for fruitful discussions 
and encouragement. 



[1] 


B. 


Grammaticos, A. Ramani, and V. Papageorgiou, Phys. Rev. Lett. 67, 1825 (1991). 


[2] 


J. 


Hietarinta and C. Viallet, Phys. Rev. Lett. 81, 325 (1998). 


[3] 


H. 


Sakai, Commun. Math. Phys. 220, 165 (2001). 


[4] 


T. 


Takenawa, J. Phys. A. 34, L95 (2001). 


[5] 


A. 


Ramani, B. Grammaticos, and Y. Ohta, Commun. Math. Phys. 217, 315 (2001). 


[6] 


A. 


Ramani, B. Grammaticos, and Y. Ohta, J. Phys. A. 34, 2505 (2001). 


[7] 


T. 


Tsuda, J. Phys. A. 37, 2721 (2004). 


[8] 


K. 


Kajiwara, T. Masuda, M. Noumi, Y. Ohta, and Y. Yamada, J. Phys. A. 36, L263 (2003) 


[9] 


G. 


Quispel, J. Robert, and C. Thompson, Phys. Lett. A. 126, 419 (1988). 


[10] 


G. 


Quispel, J. Robert, and C. Thompson, Physica. D 34, 183 (1989). 


[11] 


A. 


Ramani, B. Grammaticos, and J. Hietarinta, Phys. Rev. Lett. 67, 1829 (1991). 


[12] 


R. 


Hirota, K. Kimura, and H. Yahagi, J. Phys. A. 34, 10377 (2001). 


[13] 


T. 


Tokihiro, D. Takahashi, J. Matsukidaira, and J. Satsuma, Phys. Rev. Lett. 76, 3247 (1996) 


[14] 


1. 


Matsukidaira, J. Satsuma, D. Takahashi, T. Tokihiro, and M. Torii, Phys. Lett. A. 255 



287 (1997). 



15 



